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Abstract. We calculate and provide analytic fits of the factors which describe the reduction of the neutrino 
emissivity of modified Urea and nucleon-nucleon bremsstrahlung processes by superfluidity of neutrons and protons 
in neutron-star cores. We consider 1 So pairing of protons and either 1 So or 3 P2 pairing of neutrons. We analyze 
two types of 3 P2 pairing: the familiar pairing with zero projection of the total angular momentum of neutron 
pairs onto quantization axis, m.j = 0; and the pairing with \mj\ = 2 which leads to the gap with nodes at the 
neutron Fermi surface. Combining the new data with those available in the literature we fully describe neutrino 
emission by nucleons from neutron star cores to be used in simulations of cooling of superfluid neutron stars. 

Key words. Stars: neutron - dense matter Neutrino Emission from Superfluid Neutron-Star Cores: 
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1. Introduction 



It is well-known (e.g., Yakovlev et al. 1999| , 2001) that 
cooling of neutron stars depends on the properties of mat- 
ter in the neutron star cores. In spite of great progress in 
understanding the neutron-star physics, many properties 
of this matter are still known with large uncertainty. A 
comparison of the theoretical cooling models with obser- 
vations of thermal emission from isolated neutron stars 
gives a potentially powerful method to explore the inter- 
nal structure of neutron stars. For a successful modeling of 
the cooling one needs reliable values of neutrino emissivity 
in different neutrino reactions. 

In this paper we consider the matter of neutron star 
cores (at densities p *S 1.5 x 10 14 g cm -3 ) composed of 
neutrons (n), protons (p), and electrons (e). It is gener- 
ally agreed the neutrons and protons can be in superfluid 
state (as reviewed, e.g., by Lombardo & Schulze, 2001). 
Superfluidity affects the neutrino emission and thus the 
cooling of neutron stars. According to numerous micro- 
scopic calculations, the proton pairing occurs in the sin- 
glct ( 1 Sp) state of proton pairs. Following Yakovlev et al. 
( 1999|) we will call this pairing as pairing A. The neutron 
pairing occurs either in the state or in the triplet state 
( 3 P 2 ). Neutron pairing A takes place in the matter of sub- 
nuclear density (p po, where po = 2.8 x 10 14 g cm -3 is 
the saturated nuclear matter density), while the 3 P2 pair- 
ing is efficient at higher p. We consider the 3 P2 pairing of 



two types denoted as B and C. Pairing B occurs in a state 
of a neutron pair with zero projection of the total angular 
momentum on the quantization axis, mj = 0. This pair- 
ing has been studied in the majority of papers devoted to 
3 P2 pairing of neutrons. Pairing C occurs in a state with 
\mj\ — 2. It has been the subj ect of some studies (as re- 
viewed, e.g., by Yakovlev et al. 1999). The actual type of 



neutron pairing (A, B, or C) corresponds to the state with 
minimum free energy. Pairing C seems to be less realistic 
than B but cannot be completely ruled out by contem- 
porary microscopic theories. For example, Muzikar et al. 
( 1980| ) showed that it realizes in matter with strong mag- 
netic field (B > 10 16 G). Amundsen & 0stgaard fll985| ) 
found that the energetically preferable state of the pair 
can be a superposition of states with different mj. The 
specific feature of pairing C is that it leads to superfluid 
gaps with nodes at the neutron Fermi surface producing 
qualitatively different effect on neutrino processes than 
pairing B (or A). 



Note that we do not consider another case: 3 P2 neutron 
pairing with \mj\ = 1. In this case, just as in cases A and 
B, the superfluid gap does not have any nodes at the Fermi 
surface (e.g., Amundsen & 0stgaard 1985). Therefore, we 
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expect that the results will be similar to those for pairing B 
or A. On the other hand, the consideration of the \mj\ = 1 
pairing is technically much more complicated since the 
superfluid gap depends not only on the polar angle $ of 
neutron momentum at the Fermi sphere (see below) but 
also on the azimuthal angle ip. 
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Let us remind five main neutrino generation mecha- 
nisms in the neutron-star cores. 

(1) Direct Urea process is the most powerful neutrino 
process. It consists of two successive reactions 



(4) Neutrino emission due to Cooper pairing of nucle- 
ons (N= n or p) actually consists of neutrino-pair (any 
flavor) emission 



N — > N + v + is 



P 



P 



(1) 



(7) 



where v c and v e are electron neutrino and antincutrino. It 



is allowed (Lattimer et al. 1991) only in matter of suffi 



by a nucleon whose dispersion relation contains an energy 
gap. The process was proposed by Flowers et al. (1976) 



ciently high density (typically a few times of po) for model 
equations of state with high symmetry energy (rather high 
fraction of protons) . The reduction of direct Urea process 
by proton superfluidity A and neutron superfluidity (A, 
B, or C) was analyzed by Levenfish & Yakovlev (1994). 

(2) Modified Urea process consists of two branches. 
Two successive reactions (direct and inverse) 

n + n— ►p + n + e + t'ej p + n + e— >n + n + z/ e (2) 
form the neutron branch. Two similar reactions 



for neutron superfluidity of type A. The extension to the 
neu tron superfluidity B and C was done by Yakovlev et 
al. ( 1999 ). The case of p roton super fluidity A is described, 
e.g., by Yakovlev et al. (|l999| , |200l|) . In the absence of su- 
perfluidity, the reaction is forbidden by energy-momentum 
conservation. 

(5) Neutrino-pair bremsstrahlung at electron-electron 
scattering (Kaminker & Haensel 1999), 



n + p — > p + p 



p + p + e^n + p + z/o 



(3) 



form the proton branch of the process. The process is 
the most powerful neutrino mechanism in non-superfluid 
neutron-star cores where the direct Urea process is for- 
bidden. It (or at least its neutron branch) is open in the 
entire stellar core. 

The neutrino emissivity of this process in non- 
superfluid matter was considered by a number of authors 
(references can be found in Yakov lev et al. 1999 ), par- 
ticularly, by Bahcall & Wolf ( |l965[), Fri man & Maxwell 
1979D , and Yakovlev & Levenfish ( |l995|) . The latter au- 



thors studied the reduction of the process either by pro- 
ton superfluidity A, or by ne utron superfluidity (A or B). 
Levenfish & Yakovlev ( |l996p suggested a simple approx- 
imate method to account for the combined effect of the 
neutron and proton superfluidities. It is based on the sim- 
ilarity relations of the factors which describe the super- 
fluid reduction of the direct and modified Urea processes. 
These results were used in simulations of the neutron star 
cooling (as reviewed by Yakovlev et al. 



e + e— s-e + e + ^-l-^, (8) 

is much weaker than other processes in non-superfluid 
matter. However, it is almost independent of superfluidity 
and may be the leading mechanism in superfluid matter. 

The present paper is organized as follows. In Sect. 2 we 
present general equations for modified Urea process and 
analyze the reduction factors. In Sect. 3 we consider the 
reduction factors of nucleon-nucleon bremsstrahlung pro- 
cesses. In Sect. 4 we study the efficiency of various neu- 
trino processes in the cores of neutron stars for different 
superfluidity types. Analytic fits of the reduction factors 
of the modified Urea process are given in Appendix. 

2. Modified Urea process 

2.1. General equations 



As discussed, e.g., by Bahcall & Wolf (1965) and Friman 



& Maxwell (1979), the general expression for the neutrino 
emissivity of modified Urea process can be written as (h = 
c = ks = 1): 



1999| , |200ll ). We q = 2 / 

the reduction of J 



present a more accurate calculation of the reduction of 
the modified Urea process by combined action of proton 
superfluidity A and neutron superfluidity (A, B, or C). 

(3) The neutrino-pair bremsstrahlung at nucleon- 
nucleon scattering can be of three types: 



d 3 Pj 



d J p c 



d a p v 



2e c (2^) 3 26,(2^)3 " 



C 



(2n)U(E f -E l )6(P f ~P i )- J2\ M 



(9) 



spins 



n + n — > n + n + v + v. 
n + p — » n + p + v + v, 
p + p — > p + p + v + v, 



(4) 
(5) 
(6) 



Here, v and v stand for neutrinos of any flavor. In a nor- 
mal (non-superfluid) matter the bremsstrahlung processes 
are weaker than the modified Urea process (e.g., Yakovlev 
et al. 1999). However they may be more important in su- 
perfluid matter. If proton superfluidity is of type A, and 
neutron superfluidity is of type A or B then superfluid 
reduction of the processes can be descr ibed by the formu- 
lae presented by Yakovlev et al. ( 1999 ). In this paper we 
develop analogous description for neutron superfluidity C. 



where Pj is a nucleon momentum (J = 1, 2, 3, 4); p e and 
e e are, respectively, the momentum and energy of an elec- 
tron; and p v and e v are the momentum and energy of 
a neutrino. The delta function S(Ef — Ej) describes en- 
ergy conservation, while S(P / — P{) describes momentum 
conservation. The indices i and / refer to the initial and 
final particle states. Furthermore, C means the product of 
Fermi-Dirac functions or corresponding blocking factors of 
the nucleons and the electron; \M\ 2 is the squared matrix 
element. Summation is carried over all particle spins. The 
factor 2 in the denominator before the summation sign is 
introduced to avoid double counting of the same reactions 
involving identical particles. The overall factor 2 doubles 
the emissivity of elementary (direct or inverse) reaction 
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of the process assuming beta-equilibrium. Since neutrons, 
protons and electrons in neutron-star cores are strongly 
degenerate, one can use the phase-space decomposition 
(e.g., Friman & Maxwell 1979) which yields: 



A = An 



1=1 Ax v x v 



ri / 

5 

n 



ME* > 



(10) 

(11) 



\ j=1 



6 [ ^Xj-x v ) .(12) 



The factor A contains integrations over orientation of par- 
ticle momenta (j = 5 refers to an electron); dflj is a solid 
angle element in the direction of p^. All lengths of the 
momenta Pj with j < 5 can be set equal to the appro- 
priate Fermi momenta pp. A typical neutrino momentum 
p v is determined by the temperature T, p v ~ T <C pf- 
Thus, we neglect p v in the momentum-conserving delta 
function and integrate over orientations of p v in A (which 
gives a factor of 4m). The factor / given by Eq. ( p"2| ) con- 
tains the integrals over the dimensionless neutrino energy 
X u = Pv/T — e v /T and the dimensionless energies of other 



particles Xj 



(p - pfj)/ t ; fj = I ex p(^-) + i]" 



Finally, Eq. (|To|) contains the products of the densities 
of state of the particle species 1 < j < 5, m* being the 
effective mass at the Fermi surface. For non-relativistic nu- 
cleons (which we consider here), m* is mainly determined 
by in- medium effects. For the electrons (j=5), m* = fj,, 
where ^ is the electron chemical potential. In the absence 
of superfluidity Eq. © gives I = I a = 11513 tt 8 /120960. 
For the neutron branch ^ of th e process, Eq. ( |ll| ) yields 
(e.g., Shapiro & Teukolsky [L983D: 



2ir (Air) 4 
PL 



(13) 



This result is valid as long as p^ > PFp+P Fc - Otherwise Aq 
is given by Eq. (13) in Yakovlev & Levenfish ( |1995| ) but in 
that case the modified Urea process becomes insignificant 
because the direct Urea process dominates. 

For the proton branch (||) at p Fn > 3p Fp — p Fc one has 



while at p^ < 3p^ - p Fc (and 3p^ > p Fc ) 



(14) 



A p 




(15) 



where Q = 1 if the proton branch is allowed by momentum 
conservation (p F < 3p Fp + p Fc ), and 8 = otherwise. 
Notice that Eq. ( |15| ) can be useful if dense matter contains 
other particles but n, p, and e (e.g., muons). 



The difference of Eqs. @ and (|T|) or @ is the con- 
sequence of the fact that p Fn is significantly larger than 
p Fp in neutron star matter. 

Combining these results one can obtain the neutrino 
emissivities Qq and Qq in nonsuperfluid matter. The emis- 
sivity Qq was calculated by Friman & Maxwell (1979), 
using the one-pion-exchange approximation for calculat- 
ing the matrix element, |M| 2 , and Qq was calculated by 



Yakovlev & Levenfish (1995) using the same technique. 

Now consider the modified Urea process in the pres- 
ence of superfluidity of neutrons and protons. A onset of 
superfluidity leads to the appearance of an energy gap S 
in the momentum dependence of the particle energy e(p). 
Near the Fermi surface (\p — pp\ <C Pf), this dependence 
can be written as (e.g., Lifshitz & Pitaevskii 



£ = M - V$ 2 + V 2 P<P F > 



e = [i + \/ 5 2 + rj 2 p > p F . 



(16) 



Here, rj = v F (p — p F ), v F and p F are the nucleon Fermi 
velocity and Fermi momentum, and fi is their chemi- 
cal potential. For the conditions of our interest, 6 -C H- 
Furthermore, <5 2 = A 2 (T)F(tf), where A(T) is the gap 
amplitude, which determines the temperature dependence 
of the gap width, while is the factor which depends 

on the angle $ between the quantization axis and the par- 
ticle momentum. The functions A(T) and F($) depend 
on superfluidity type (e.g. Yakovlev et al. 1999[ ). For cases 
A, B and C: 



F A (0) = 1, T cA = 0.5669 A(0); 

F B (??) = 1 + 3cos 2 tf, T c b = 0.8416 A(0); 



F c (tf) = sin 2 ^, 



T cC = 0.4926 A(0). 



(17) 
(18) 
(19) 



The gap amplitude A(T) is assumed to be governed by 
the standard equations of the BCS theory (e.g., Tamagaki 
1970D ; A(0) is related to the critical temperature T c as 



indicated above. 

For further analysis it is convenient to introduce the 
dimensionless variables: 



£ — /X 

T 



— sign(x)\Jx 2 + y 2 , 



1 

T' 



y = 



A(T) 
T 



(20) 



While T decreases from T c to 0, the parameter v varies 
fro m to oo as described, e.g., by Eq. (11) in Yakovlev et 
al. (tL999|) . 

We assume that the neutrino emissivity in superfluid 
matter can be calculated from Eqs. (|Io| ) — ( |l2|) by replac- 
ing Xj — > Zj for all particle species which are in superfluid 
state. This assumption is widely used in the literature; its 
validity is discussed by Yakovlev et al. (|2001 ). In this ap- 
proximation, the neutrino emissivity of the modified Urea 
process can be written as 



Q n = QoR a , Q p = Q p R p , 



(21) 



where Qq and Qq are the emissivities in a non-superfluid 
matter, while R n and R p are the factors which describe 
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the reduction of the emissivity by nucleon superfluidity 
(i? N < 1). Generally, these factors can be written as 



Jn 

rN AN ' 

1 a A o 



/ 5 n CO 

Y[dSlj / dx v xl 




dxj f(zj) 



(22) 



The notations are the same as in Eqs. jlO| ) — (|l2"|). 

We have composed a code which calculates the reduc- 
tion factor ( p2|) for proton superfluidity A and neutron 
superfluidity A, B, or C. The code has been tested by 
comparing with the analytical asymptotes at large v\ and 
v 2 and with the results of Yakovlev & Levenfish (1995) 
who considered superfluidity of either protons or neutrons. 
The results have also been compared with those calculated 
from Eq. (H|) under simplified assumption p^, p Fc <g; p a 
discussed below (see Eq. [y]). 

Notice that the results of this section can also be used 
to describe modified Urea pr ocess with muons instead of 
electrons (see Yakovlev et al. 2001, for details). 




Fig. 1. Four regions of v\ — v n and v 2 = v p where the 
reduction factors i? N of the neutron and proton branches 
of modified Urea process can be fitted by different expres- 
sions. In regions I, II, and III at v\ 3> 1 and v 2 3> 1 the 
factors R n of neutron branch have different asymptotes. 



2.2. Reduction by superfluidity of neutrons and 
protons of type A 



In this case Eq. (|22|) can be simplified. For pairing A, the 
dimcnsionless energy gap y A is angle- independent. This 
allows one to decompose the integrals over the angles and 
over the dimensionless energies Xj . For the neutron branch 
of the modified Urea process we get 



#AA 



G(a) = 



120960 



115137T 8 



+ 00 



n / a ^ 

x G(z\ + z 2 + z 3 + z A ), 
{x + af 



-t-oo 



dx 



1 



(23) 
(24) 



where j = 1,2, and 3 enumerates the reacting neutrons, 
while j = 4 refers to a proton. 

Let here and hereafter v\ = v n refer to neutrons, and 
V2 = v p refer to protons. The reduction of the proton 
branch is evidently given by 

R p AA (vi,v 2 ) = Rl A (v 2 , Vl ) (25) 



It is useful to find the asymptotes of R AA from Eq. (|23|) 
for strong superfluidity i.e., for large values of vi and v 2 - In 
this case we can introduce three regions of the parameters 
(Fig. 1), where the asymptote of R AA has different forms. 
Region I corresponds to V\ > v 2 ; region II corresponds 
to v 2 > vi > v 2 /3; and region III corresponds to v\ < 
v 2 /3. The asymptotes presented below are valid at (vi — 
v 2 ) 2 3> V\ and (v 2 — 3vi) 2 ^> v 2 , i.e., not too close to the 



boundaries between regions I and II and between regions 
II and III. 

For example, we outline the derivation of the asymp- 
tote of R AA from Eq. (|23| ) in region I; the derivation in 
other regions is similar. Clearly, the integral ( ^3|) can be 
subdivided into several parts in such a way that any sin- 
gle part contains integrations from — oo to and/or from 
to +oo. Now let us introduce the convenient notations 
for these parts. Let R(2, — 1) mean a five-dimensional in- 
tegral containing the integration from to +oo over two 
neutron variables, and over — oo to over a proton vari- 
able (in this case, the integration over the third neutron 
variable is assumed to extend from — oo to 0). Splitting 
the initial integral ( p3| ) into the elementary integrals,we 
see that the same integral R(2, — 1) enters the sum three 
times. Thus, it is sufficient to calculate the integral once 
and multiply by 3. 

In this way we obtain eight integrals of differ- 
ent types: i?(3,+l), #(3,-1), R(2,+l), i?(2,-l), 
R(l,+1), i?(l,-l), R(0,+1), and #(0,-1). In the 
limit of strong superfluidity (vi t2 ^> 1), each of 
them is exponentially small. The exponentials are: 

i?(3, +1) oc exp(— 3wi — v 2 ), R(l, +1) oc exp(— 2vi), 
R(3, —1) oc exp(— 3ui), R(l, — 1) oc exp(— 2«i — v 2 ), 
i?(2, +1) oc exp(-2«i - v 2 ), R{0, +1) oc exp(-3ui), 
R(2, -1) oc exp(-2ui), #(0,-1) oc exp(-3«i - v 2 ). 

It is seen that the main contribution into the asymp- 
tote comes from R(2, — 1) and i? (1,4-1). These terms 
have the same exponential but i?(2, —1) has a larger pre- 
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exponent. Therefore, it is R(2,—l) which gives the main 
contribution in region I: 



#(2,-1) = 



dxi f{zi) dx 2 f(z 2 ) 



120960 

115137T 8 

•0 

x | / dx 3 f(z 3 )dx 4 f(z 4 )G(a), (26) 



o 



where a = zi+z 2 +z 3 +z 4 . One has G(a) — ► a 4 /4 as a — > 00 
and G(a) — ► 6exp(a) asa~> —00. We are especially inter- 
ested in large and positive a for which G(a) » a 4 0(a)/4, 
where 0(a;) is the step function. In this approximation, 
Eq.(p6|) can be rewritten as: 

120Q60 f f + °° 
#(2,-1) = JYbm^E J I dx i f(z 1 )dx 2 f(z 2 ) 



dx 3 f(z 3 )dx 4 f(z 4 )—9(a). (27) 



Two integrals in Eq. ( p7| ) over the neutron variables 
X\ and X2 are rapidly converging because of the presence 
of exponentially decreasing functions. Accordingly, the re- 
gion of space, where these two variables produce the main 
contribution, is rather small. Thus, it is sufficient to set 
Z\ = v\ and z 2 = v\ in the function. Then the inte- 
grations over x\ and x 2 are separated and done analyti- 
cally. Neglecting exponentially small terms while integrat- 
ing over X3 and x 4l we obtain the asymptote of R\ A in 
region I: 



#aa — 3 



120960 7T 



1 V\ 



h(a) 



where r 



11513tt 8 2 

i/r 

x 0(2 - r 3 - ar 4 ), < a < 1 



da; 3 da; 4 (2 — r 3 — ar 4 ) 4 



(28) 



(29) 



1. 



In the same manner in region II we obtain: 



^AA 



h{a) 



120960 tt 13/2 i/ 2 _ V1 _ V2 

— ' ''' '' L'l 



11513tt 8 2 1 2 



dxi dx 2 (1 + a — rx— r 2 ) 4 



115137T 8 

i/r 

X 0(l + a-n -r 2 ), 1 < a < 3, 



(30) 



(31) 



if = 



V9 f 2Q (486^ + 747, 2 , 2 + 16, 4 ) 



3 2 /o 2 , or 2. , 3t>i + y/9« 2 - v\ 

-—vi v 2 (3v 2 + 3bv 1 ) In . (32) 



''2 



In region III 

120960 /7T\i/2 



^ A ~ 115137T 8 



/7T\ 1/2 3 9/2 -v 2 j ( V l\ 



(33) 



/3(a) 



+ OC 



dx\ dx2 dx 3 



4 j j .1,, 
x 0(a - ri - r 2 - r 3 ), 



1 (ri + r 2 + r 3 ) 



(34) 



where 3 < a < 00. 

The asymptotes themselves contain complicated inte- 
grals. Thus we have calculated the asymptotes of the in- 
tegrals (E9f) and (p54). In region I we have: 



0.027570965 



h = 



h = 0.0785398 (1 - af 



0, 
1. 



(35) 



The asymptote of I 2 {a) at a — > 1 can be determined by 
matching the reduction factors at the boundary of regions 
I and II. 

In region III we have 



I 3 = 0.214307 

1 a 



1-1 

a 



11/2 



3. 



840 



a 



00. 



(36) 



We have numerically calculated the functions Ii(a), 
/2(a), and /3(a) and proposed analytic fits which repro- 
duce the results and the asymptotes (p9|), (fTlj), and (34) 
with the maximum error less than 1%. For region I, the 
fit is: 



h 



0.02757096 (1-a? 1 ) 5 



a (1 + 30a 2 



60a 4 + 107.186a 6 )^ ' 
0.1684. 



where p\ = 1.473 and p 2 
For region II: 

h pa 0.05 (a pi - l) 5 (a P2 -p 3 ) 2 , 

with p! = 1.210, p 2 = 0.222, and p 3 = 0.215. 
For region III: 

5.5 

h - — I' - — I (I 



a. 

840 



i-4 

or 



Pi N 



(37) 



(38) 



(39) 



where pi = 8.363, p 2 = 1.427, and p 3 = 1.8978. 
The fits of i?AA are gi ven m Appendix. 

2.3. The neutron branch reduced by proton 
superfluidity A and neutron superfluidity B 

Now consider neutron superfluidity of type B. According 
to Eqs. ([20] ) and (|l8j), the dimensionless energy gap of the 
neutrons, ye, is angle-dependent. The asymptotes of R\ B 
in this case can be obtained from Eq. (|22]). As before, j=l, 
2, and 3 enumerates neutrons while j—A refers to a proton. 
Equation (p2h can be written as 



1 



PL 



120960 

11513tt 8 (4tt) 3 p Fc p Fp 

X /6l (Pphl+Pph2+PPh3 



df7idfi 2 df7 3 



IPFnl + PFn 2 + PHiJ 



n 



+ OC 



da;,- /(*,•) 



G(zi + z 2 + z 3 + z 4 ), (40) 



where G(a) is given by Eq. (|24|); is the step function: 

e = 1 if b Fc -PFpl < IPf„i +Pf„ 2 +P&3I ^ Pfc +Pf p - and 
= otherwise. 
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Calculations show that the reduction factor is al- 
most insensitive to variations of particle Fermi-momenta. 
Let us obtain a simplified expression for i?^B by setting 
p Fp =p Fc =0 in Eq. ( p^ ) and integrating over dO p and dfl e . 
The result is 



A 



AB 



PL 



120960 



11513tt 8 2(2tt) 5 



3 

E 



pj 



-\-oo 

n / d ^ fa) 

i=i 



G(zi + z 2 + z 3 + z 4 ). (41) 



As for neutron superfluidity A, we have different asymp- 
totes of R\ B in regions I, II, and III. For example, con- 
sider the asymptote at large v\ and v 2 in region I. It 
is easy to see that the main contribution comes from 
#(2,-1) (see Eq. ©). In Eq. @ it is sufficient to 

2, 3) and 



sign(xj)^/^ +vf(l + 3c 2 ) (for j 



=1, 



SCt Zj = 

Z4 = sign(x4) A /a?4 + ijf. Here, Cj = cosi?j. The simplified 
reduction factor will then be rewritten as: 



R1b " 2(2tt) 2 



J dc 1 dc 2 dc 3 D(c 1 ,c 2 ,c 3 )R(2, 
D(ci,c 2 ,c 3 ) = / dipid(p 2 dtp3 8 V" — 

^0 Vfe^y 



-1), (42) 



4tt 6>(0.75 - cic 2 - cj - cj) 



ClC 2 



5(ci + c 2 + c 3 ). 



'1 L 2 

Notice that the main contribution into ( [42| ) comes from 
the range of angles ci « c 2 ~ 0. Then it is sufficient to set 
Zi w v\ + 0.5a;f 1 + 1.5c 2 vi (i = 1, 2) in the exponentials 
in R(2, —1), and z\ = z 2 = «i in all other functions. This 
leads to the following asymptote in region I: 

120960 n , -** Ix (?A t (43 ) 



#A 



AB 



11513tt 8 3%/3 



where I\(a) is defined by Eq. (29). Now let us intro- 
duce three functions J,-(ci, c 2 , C3, a) (i = 1, 2, and 3) , 
which formally coincide with those given by Eqs. fl29J), 
(^l|), and (HJ) with the only difference that now rj — 

1, 2, and 3. 



3c 2 , j-- 



Then, the asymptote in region II will be written as 



#AB 



120960 n- 1 / 2 



6„ 1/2 



11513tt 8 8x/6 
< y dcidc 2 L>(ci,c 2 ,0)/ 2 ^ci,c 2 ,0, 



120960 7T 3 / 2 



,1/2 



-3vi 



11513tt 8 3\/6 
where K is given by Eq. (|3 
120960 71-- 3 / 2 



A'. 



(44) 



In region III we obtain: 



R\ 



AB 



115137T 8 



8x/2 



3 9/2 -Do 

yi« 2 e 2 



x / dcidc 2 dc 3 D(ci,c 2 ,C3)/ 3 ci,c 2 ,c 3 . 



''2 



,(45) 



We have numerically integrated R\ B from Eq. ( flGf ) for 
a dense grid of «i and w 2 . The calculations have been con- 
ducted at p Fp = 0.11 p Fn and p Fc — 0.1 p Fn . As mentioned 
above, the reduction factor is rather insensitive to varia- 
tions of these parameters. The variations of i?^ B to the 
changes of the particle Fermi momenta within reasonable 
limits (pfc, p < (0.3 - 0.4)p Fn ) obtained in some test runs 
are of the order of estimated error of numerical integra- 
tion. The fits of R\ B are given in Appendix. 

2.4. The proton branch reduced by proton superfluidity 
A and neutron superfluidity B 

In this case Eq. (p2|) can be simplified as 



•f. A D 



120960 1 



AB 11513tt 8 2 J_ x 



dc 4 7(zi,z 2 ,z 3 ,z 4 ). 



(46) 



Here, C4 = cos $4. The function I(zi, z 2 , z 3 , Z4) is defined 
by Eq. (|4(]) with the only difference that now j = 1,2, and 
3 refer to protons, while j ' = 4 refers to a neutron. 

One can easily obtain the asymptotes of the reduction 
factor at large values of v\ and v 2 . For the proton branch of 
modified Urea process, the regions where the asymptotes 
are different can be found from neutron-branch regions by 
replacing v\ =^ v 2 . For instance, at v 2 > V\\ 



AB 11513tt 8 4 2 



J dc 4 l x ( c ^) » ( 47 ) 



where ii(c4,a) is defined by Eq. ( |29| ) with r 4 
sjx\ + 1 + 3cf . At vi>v 2 > Vi/3: 



7? p 



120960 /7TN3/2 



AB 11513tt 8 V2 



13/2 - Vl -i 



i'2 



120960 ^ 2 n /2 r 1 
115137T 8 29/ 2 WlW2 6 ' <l " 



dx4 (3 -r 4 ) ( 3 -f4 

V "2 / V V 2 



At w 2 < vx/3: 

_ 120960 7T 
AB ~ 11513tt 8 27! 



(48) 



(49) 



The fits of -R^b are gi ven i n Appendix. 

2.5. The neutron branch reduced by proton 
superfluidity A and neutron superfluidity C 

The most important feature of this case is that the energy 
gap vanishes at the poles of the Fermi sphere (see Eqs. ( 19 ) 
and (pp[)). Equation ( p7j| ) remains valid in this case. The 
calculations of R\ c have been done at p Fp = 0.11 p Fn and 
p Fe = 0.1 p Fn . As in the previous cases, the reduction fac- 
tor is rather insensitive to variations of these parameters. 
The results are approximated by the expressions given in 
Appendix. 
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2.6. The proton branch reduced by proton superfluidity 
A and neutron superfluidity C 

Equation ( |46| ) remains valid in this case. Since j/4 = yc = 
V\ sin $4 vanishes at the poles of the Fermi sphere, the re- 
duction factor R AC varies with v\ as a power-law (rather 
than exponentially). It is easy to determine its behavior 
at large V\. One can see that in this case the main con- 
tribution into integral (|46| ) comes from the region where 
sin 1^4 <C 1. Thus, we have 

120960 2 f +OC 

RAC = 115137r 8 V l J o dt 1 T ( Zl ' Z2 > Z3 ' Z4 )' ( 50 ) 

Here, I is given by Eq. (^), in which Z4 = 
sign (X4 )y/ x\ + 1 2 . The difference of exact Eq. ( |46| ) from 
Eq. ( p0| ) is less than 2% at v\ 25. As in the previous 
cases, R\q has been calculated numerically. For vi > 25, 
we have used Eq. ( |50| ) and fitted our results by a simple 
equation 

dp -2 f Pi u 2 \ 

pa = 0.289297, p 2 = 0.152895, p 3 = 0.228536, 

Pi = 2.250496, p 5 = 3.885153- 10~ 3 . (51) 

The fits of -R^c f° r Vl — ^5 are gi ven m Appendix. 




2 4 6 8 10 



v 

Fig. 2. Reduction of the neutrino emissivity by neu- 
tron and proton superfluidities of type A in the neutron 
branch of the modified Urea process versus v (v% — v sin tp, 
V2 = vcosp) at ip — 0, 45, 63, and 90°. Solid lines show 
our results, and dashed lines are ob tained from similarity 
criteria, e.g., Yakovlev et al. ( |l999| ). 



3. Approximate reduction factors of the 
NN-bremsstrahlung in superfluid matter 

Now consider the superfluid suppression of the neutrino- 
pair emission in the nucleon-nucleon bremsstrahlung pro- 
cesses (||) - (H). In the absence of superfluidity the 
emissivities of NN-bremsstrahlung processes in the one- 
pion-exchange approximation are given, for instance, by 
Yakovlev et al. ( |l999[ ). 



In analogy to Eq. 
reduction factors R N 

Qnn = <2nno-R NN > 
where 



rNN = 945 P Hl 



Ij), one can introduce the superfluid 



(52) 



NN 



rNN 



164tt 8 32 



n 

i=i 



dxj f(zj) 



X S(zi + Z 2 + Z3 + Zi ~ x v ). 



(53) 



Here, Zj is given by Eq. (p0[), and j enumerates initial- 
state and final-state particles participating in the re- 
actions. Accordingly, if we consider neutron superflu- 
idity C, and proton superfluidity A, then we have 



Zj = sign(xj ) ^Jx 2 + v\ sin 2 
sign(xj )\Jx 2 +V 2 for protons. 



j for neutrons, and Zj 



The factor R P A P X was accurately calculated by Yakovlev 
& Levenfish ( 1995 ). For the neutron- proton process we 
suggest the similarity relation of the form 



DIP 

^AC 



KIM 



(54) 



Here, R^ c (vx, V2) is the reduction factor of the direct Urea 
process determined by Levenfish & Yakovlev (1994). The 
subscript pA means that superfluidity of protons is of type 
A (and neutrons are non-superfluid) . The factor R^ A was 
calculated by Yakovlev & Levenfish ( |1995 ). 

An analysis of i?"^ for neutron-neutron brems- 
strahlung, Eq. (Q), is more sophisticated (since no simi- 
larity criterion can be formulated). Let us study the re- 
duction factor R^q at large Vi from Eq. (^3|). Now j = 1—4 
refer to neutrons. One can see that the main contribution 
to comes from the range of angles sin-dj ss 0. Since 
the sum of the Fermi momenta of reacting neutrons must 
be equal to zero, the Fermi momenta should concentrate 
to the poles of the Fermi sphere: two momenta to one pole 
and other two momenta to the other pole. Now we expand 
all functions in series over 0j and integrate over tp4, <p 3 , tpi, 
and $4. In this way we obtain 



945 96tt 
nC ~ 164^8 32^3 



d0id0 2 d0 3 ft P 



x 0(0* + 0^ - 0? > 0) 



dt 



^{\-t 2 ){m 2 -t 2 ) 



0(t < m), 



(55) 
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#1 



02 03 



(56) 



In Eq. (|55|) we introduce t = cos ip 2 ; I nn is the same 
as in Eg. (|5^) with the only difference that now $4 = 



— -# 2 . Introducing yj — Vi'&j, j — 1,2,3, and 



taking into account that v\ is large, we can extend the 
upper limit of integration over yj to infinity. As a result, 
the asymptote of becomes 

945 96tt r , , „ 

dyidy 2 dy 3 y 1 I 



pnn 



164tt 8 32tt 3 1 

Kvl + vi ~ Vi > 0) 
dt 



= 6(t < to) 

/o y/(l-t*)(m? -i 2 ) " 

Evaluating this integral, we obtain: 

11.533 



l iiC 



(57) 



(58) 



Let us derive an approximate formula for R^c at in- 
termediate values of v\. For this purpose we substitute 
| (yi + V2 + 2/3 + Va), where y 3 = V\ sintf,, in the argu- 
ment of the function R^\(vi) 



(as described by Yakovlev 



et al. 1999) and integrate this function over $ 



nun 



n 



dsini?,- 



#Sa( J(l/i+«2 + »a + y4)),(59) 



where Pl = 0.442995, p 2 = 0.250953, p 3 = 0.410517, p 4 = 
7.09171T0 -3 . The calculation and fit error does not exceed 
2%. 



4. Discussion 

The results of Sect. 2 allow us to compare the exact and 
approximate reduction factors of modified Urea process. 
The comparison is illustrated in Figs. 2, 3, and 4. The 
figures show the dependence of the calculated reduction 



factors, i?AA> ^AB' an d ^ab> on v 



where / = 2.248 is chosen to satisfy the asymptote (p8|). 
For vi < 25, has been calculated numerically. At 

v± J> 25 our numerical results agree with the asymptote 
( |58| ) within 1%. The numerical results for V\ < 25 can be 
fitted as 



i nC w exp 



(1 



P2 v f 



(60) 



«2 at sev- 



eral values of ip (ip is the polar angle in the vi-v 2 plane; 
tan(tp) = (v\/v 2 )). Our results (solid lines) are compared 
with the approximate reductio n fact ors (dashed lines) con- 
structed (e.g., Yakovlev et al. 1999] ) using the criteria of 
similarity between the reduction factors for different neu- 
trino reactions. The approximate factors have been used 
in a number of simulations of neutron star cooling. One 
can see that the difference of the approximate reduction 
factors from the exact ones increases with increasing v 
(but for ip = 0° in Figs. 2, 3 and ip = 90° in Fig. 4). 

Now let us answer the question which neutrino gener- 
ation mechanism dominates in a superfluid neutron-star 
core. Taking into account the above results we can calcu- 
late the emissivities of all main neutrino processes (Sect. 1) 
for proton superfluidity A and any neutron superfluidity, 
A, B, or C. Figures 5 and 6 show which process dominates 
at different values of T cn and T cp . Figure 5 shows the ef- 
fect of neutron superfluidity B, while Fig. 6 - the effect 
of neutron superfluidity C. Both figures are plotted using 
an equation of state suggested by Prakash et al. ( 1988 ) 
(their model I of symmetry energy with the compression 
modulus of saturated nuclear matter K — 240 MeV). 
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erf 
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Fig. 5. Regions of T cn (of type B) and T cp (of type A) where the different neutrino reactions dominate. 



Three left panels of Fig. 5 illustrate standard neu- 
trino emission at p = 2po (direct Urea process forbid- 
den) for three values of the internal stellar temperature, 
T = 10 8 ,3 x 10 8 , and 10 9 K, while three right panels 
correspond to neutrino emission enhanced by direct Urea 
process at p = 5po for the same T. Figure 5 is alm ost the 
same as obtained earlier by Yakovlev et al. ( 1999 ) for an- 
other equation of state using the approximate reduction 
factors of modified Urea process. The selected values of 



T cover the temperature interval most important for the 
theory of neutron star cooling. The figures are almost in- 
dependent of p (and of equation of state) as long as p does 
not cross the density threshold of opening direct Urea pro- 
cess. One can see that if the neutrons are superfluid alone 
and T -C T cn , then the bremsstrahlung due to proton- 
proton scattering becomes dominant. If protons are su- 
perfluid alone and T <C T cp , then the main mechanism 
is neutrino emission in neutron-neutron bremsstrahlung. 
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Fig. 6. Regions of T cn (of type C) and T cp (of type A) where the different reactions dominate. 



Neutrino emission due to Cooper pairing of neutrons al- 
ways dominates at T ~ 0.4T cn provided direct Urea pro- 
cess is forbidden. If the direct Urea process is allowed 
then the Cooper-pairing neutrino emission may dominate 
provided nuclcons of one species are strongly superfluid 
while nucleons of the other species are moderately super- 
fluid. Finally, in the presence of strong superfluidity of 
protons and neutrons, all the processes involving nucleons 
are so strongly reduced that the neutrino-pair emission in 



electron-electron bremsstrahhmg dominates. Figure 6 dif- 
fers from Fig. 5 mainly by the increase of the efficiency of 
neutrino emission due to Cooper pairing of neutrons (at 
p = 2po) and direct Urea process (at p = 5po)- 

It is well known that only one neutrino process domi- 
nates at a given density in a non-superfluid neutron-star 
core. It is either direct Urea or modified Urea process. The 
situation is drastically different in superfluid matter. As 
seen from Figs. 5 and 6, Cooper-pairing neutrino emis- 
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sion becomes dominant in the presence of a weak neutron 
superfluidity. With the increase of T cn in the superfluid 
regime, modified Urea process becomes unimportant and 
Cooper-pairing neutrino emission dominates. 

5. Conclusions 

We have calculated the factors which describe the reduc- 
tion of the neutrino emissivity in the neutron and pro- 
ton branches of modified Urea process by superfluidities 
of neutrons and protons. We have considered singlet-state 
pairing of protons (pairing A) and either singlet-state or 
triplet-state pairing of neutrons (A, B or C). The reduc- 
tion factors are fitted by analytic expressions presented in 
Appendix to facilitate their use in computer codes. 

We have also considered the reduction of neutrino 
bremsstrahlung due to neutron-neutron and neutron- 
proton scattering by proton superfluidity A and neutron 
superfluidity C. We have constructed the approximate re- 
duction factors and fitted them by analytic expressions. 
We have determined also the dominant neutrino emission 
mechanisms in a neutron star core at different values of 
the critical temperatures of the neutron and protons, T cn 
and T cp , for the cases of neutron superfluidity of type B 
or C. 

Our results combined with those known in the litera- 



Table 1. Fit coefficients pi for the reduction factor -R^Ai 



ture (e.g., Yakovlev et al. |200lj) allow one to calculate the 
neutrino emissivity in a neutron-star core in the presence 
of proton superfluidity A and neutron superfluidity A, B, 
or C. The results can be useful to study thermal evolution 
of neutron stars, first of all, cooling of isolated neutron 
stars. Our cooling simulations based on the present re- 
sults will be published elsewhere. 
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In Fig. 1 we plot four regions, I, II, III, and IV, in the 
V\-V2 plane where the fit expressions for the neutron and 
proton branches of modified Urea process are different. 
This selection of the regions is the same for any type of 
neutron superfluidity, A, B, or C. 

We have calculated the reduction factors of the modi- 
fied Urea process from Eq. (^2|) as described in Sect. 2. 
Introducing the polar coordinates (v\ = v sin 95, = 
v cos ip), in regions I, II, and III we fit the numerical results 
by the expression 
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Av 2 



while in region IV we use the fit of the form 
i? N = Cexp(-A/B), 



(Al) 



(A2) 



In regions I, II, and III, the functions A, B, and C depend 
on ip. In region IV, they depend on v\ and V2- 

In the case of neutron and proton superfluidity A for 
neutron branch of modified Urea process we get the follow- 
ing fits. 
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In region I: 
A = pi+ p 2 t 2 

B = p 6 + Pit' + 

C = Pll 



Pi 



l+P3t 
2 + P9 

l+P&t 
Pl2 



P5t, 

Pmt, 



y (1 +P13Z 2 ) 3 ' 
t = cos 2 (ip), y = sin 2 (ip +p 15 ), 
z = cos 2 (ip + p 14 ) . 

In region II: 

2 , Pi 



(A3) 



A = pi + p 2 z 



1 +P3Z 



P5Z, 



B = pe 



P7Q 



Vpiiqt , 
Piny 3 



1+put 2 

C = pia+puy 3 + 1 """" — J -pny" 
1 +Pisz z 

z = cos 2 ip, t — cos 2 {<p+p%) 



q = sm 2 (p g ip + pio), y 
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A = Pi - 

B = Pio 
C = pie 

t = 



piy 

1 + P2y + pzy 2 

Piy 
l +psy + pgy 

p^y 
1 +puy + Pi5y 2 
, i 1ixp> 



3tt 

snr If+Y 



Pby + Peyt, 
2 rpny + pi 2 yt, 



pny + pisyt, 



0.321750554 



y 



sin 2 <p. 



(A5) 



In region IV: 

A = Piv\+p 2 vl +pzv\vl +piv\ +p 5 v 2 , 
B = l+p 6 vj+p 7 v 2 +p 8 Vi, 

C = 1+ P9 V 4 2 , 



(A6) 



Coefficients pi which enter Eqs. (|A3|)-( A6) are listed in 
Table 1. The calculation and fit errors do not exceed 10% 
for R\ A > 10" 5 . 

In the case of proton superfluidity A and neutron su- 
perfluidity B for neutron branch of modified Urea process 
we get the following fits. In region I: 



A = pi+p 2 t 2 
B = pe+ p 7 t 2 
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Pot 

l+ P8 t 2 
P\ 2 
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Table 2. Coefficients pi for R\ B 



i 


I 


II 


III 


IV 




1 


-0.719681 


-6.475443 


0.316041 


0.565001 




2 


-0.024591 


-1.186294 


-289.2964 


0.087929 




3 


0.297357 


0.591347 


2480.961 


0.006756 




4 


1.260056 


6.953996 


-268.8219 


1.667194[- 


-4] 


5 


0.100466 


3.366945 


1984.115 


3.782805[- 


-6] 


6 


0.148464 


-9.172994 


3503.094 


0.173165 




7 


0.253881 


-2.675793 


0.331551 


1.769413[- 


•5] 


8 


140.3699 


1.053679 


-0.265977 


7.710124[- 


-8] 


9 


0.132615 


10.38526 


1098.324 


0.001695 




10 


0.280765 


7.138369 


65528.01 






11 


0.375796 




0.024500 






12 


-0.096843 




0.120536 






13 


3.100942 




89.79866 






14 


0.275434 




5719.134 






15 


0.330574 




285.8473 






16 






0.402111 






17 






16657.19 







(A4) c = P6 + P7 z 2 + 



P<> 



■Pwz, 



1 +P8Z 

Z = cos 2 (p, t — cos 2 (ip + p s ) , 

q = sin 2 (p 9 (p+p w ), y = sin 2 (^ + 3ir/4). 

In region III: 

Pi 2 (1 + Pi5^ 2 + Pi6^ 3 + PnV 4 ) 



(A8) 



A 

B = pii 
C 



1 +Pl3</? 2 +P14V 3 
P7 

1 + p 8 ip 2 + P9<P 3 + PlO^ 4 ' 

Pl (1 + p 4 <^ 2 + P5<^ 3 + P6</5 4 ) 



1 +p 2 ip 2 +P3ip 3 



(A9) 



In region IV: 



A = piv 1 +p 2 v 2 +p 3 v 1 v 2 +P4«i, 
B = 
C = 1 



\P6«i 



p 7 v\ 



P5V1 



Psvf, 



Pov\- 



(A10) 



Coefficients pi which enter Eqs. ( A7)-( Alp| ) are listed in 
Table 2. The calculation and fit errors do not exceed 5- 
15% in those cases in which R\ B > 10~ 5 . 

In the case of proton superfluidity A and neutron su- 
perfluidity B for proton branch of modified Urea process 
we get the following fits. 

In region I: 



Po 



A = pi 
B 
C 

y = sin 2 (y>+pi5) 
In region II: 



P2V + 

Pl2 - Pl3^ 



p 4 ip 



(l+p 3 t<p) 2 
P9<f 



p 5 ttp 2 , 



(1 + p s tip + p w yt) 2 
Pio 



P11V 



(1+Pi 4 i 2 ) 2 

t = COS 2 ip. 



Pnttp, 



(All) 
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Table 3. Coefficients Pi for R 



AB 



Table 4. Coefficients pt for R\ c 



i 


I 


II 


III 


IV 




i 


I 


II 


III 


IV 




1 


0.288203 


0.398261 


0.387542 


0.272730 




1 


0.897393 


-3.471368 


0.322115 


0.175090 




2 


-0.124974 


-0.054952 


-195.5462 


0.165858 




2 


-0.045357 


-0.133540 


-15.05047 


0.088159 




3 


17.39273 


-0.084964 


3032.985 


0.005903 




3 


0.309724 


0.143230 


112.9733 


3.055763[- 


•3] 


1 


0.083392 


-0.036240 


-189.0452 


2.555386[- 


-5] 


4 


-0.739962 


3.634659 


-13.79012 


3.984607[- 


•7] 


5 


0.059046 


-0.168712 


3052.617 


2.593057[- 


7] 


5 


0.222597 


0.496579 


128.3156 


5.591497[- 


•8] 


6 


0.028084 


-0.704750 


442.6031 


0.023930 




6 


0.032104 


0.030609 


39.82789 


0.046496 




7 


-0.019990 


-0.066981 


0.041901 


0.006180 




7 


-0.054011 


0.005056 


0.164614 


1.452790[- 


5] 


8 


28.37210 


1.223731 


-0.022201 


1.289532[- 


-5] 


8 


61.73448 


0.438608 


49.07699 


4.505614[- 


8] 


9 


0.244471 


0.363094 


5608.168 


0.005368 




9 


0.195679 


-2.970431 


-3.145006 


1. 779724 [- 


•3] 


10 


—0 610470 


—0 357641 


-10761.76 






10 


—0 001851 


284703 

v.' ■ — v.- 1 I V 7 • J 


5132 076 


2 136809 - 


-41 


11 


0.023288 


0.869196 


0.064643 






11 


0.482581 


0.898355 


0.018737 


5.365717[- 


-41 


12 


0.475196 


-0.364248 


0.296253 






12 


-0.001637 


-0.036420 


0.100223 






13 


-0.180420 


2.668230 


106.3387 






13 


-0.685659 


0.407393 


4.055407 






11 


25.51325 


-0.765093 


-75.36126 






11 


1.528415 


-0.058942 


390.6242 






15 


0.281721 


-4.198753 


84.65801 






15 


-0.053834 


0.605413 


6.594365 








n nsnzisn 

u.uou^ou 












-0.452426 


2.851209 


175.7396 






17 


-0.191637 




-86.76801 






17 


-0.053502 


-0.800218 


441.3965 


















18 




1.497718 












P4 + 








19 




1.476375 








A = 


Pi +P2Z 2 + 


P5<P, 




















1 


















B = 


p 6 +p 7 z 2 + 


P8 + 








C = 


Puq c 


f2 








1 + p 9 z 












C = 


Pll -P12<P 2 


Pl-iV 2 








z = 


cos 2 <f, tb 


= cos 2 (<yS +p s ), 






1 +Pi4^ 






Qb = 


sin 2 (p 9 ip + 


Pw), t c = 


cos 2 (^ + pi 5 ), 





cos 2 ip. 



(A12) 



sin 2 (p 16 ip +p 17 ). 



(A15) 



In region III the factors A, B, and Chave the same form 
as in Eq. ([A|). 
In region IV: 



In region III the factors A, B, and Chave the same form 
as in Eq. (A9). In region IV: 



A = pivl + p 2 v\ + p 3 vjv 2 + p A v\ +p 5 vf, 
B = l+p 6 vl+p 7 v 2 +p 8 V2, 



A = p x v\ + p 2 vj + p 3 vjvj +P4vf, 
B -- 



C = l+p 9 V$. 



(A13) 



'1 +p 6 V 2 +p 7 vl +p 5 vf + 

C = l+p 9 vl+p w vl+pnvlvl. 



Psvf, 



(A16) 



Coefficients p, in (AS), ( A11 )-( [A13| ) are given in Table 
3. Calculation and fit errors do not exceed 8% for those 
values of v\ and v 2 for which R\ B <; 1CP 5 . 

In the case of proton superfluidity A and neutron su- 
perfluidity C for neutron branch of modified Urea process 
we get the following fits. 

In region I: 



Coefficients pi in (Al4)-(Al£) are given in Table 4. 
Calculation and fit errors do not exceed 10% for those 
values of v\ and v 2 for which i?^ c JeJ 10~ 6 . 

In the case of proton superfluidity A and neutron super- 
fluidity C for the proton branch of modified Urea process 
at v\ < 25 we get the following fits 

In region I: 



A = pi+p 2 t 2 + 
B = p 6 + P7 t 2 + 
C = p n 



Pi 



l+p 3 t 

P9t 

(l+Pstq) 2 
P12 



Pit, 



Piot, 



P17S 2 , 



s(i+ Pl3 z^ 

t = COS 2 ip, Z — COS 2 (<y9 + P14), 

q = sin 2 (^+pi 6 ), s = sm 2 ((p+p 15 ). 
In region II: 
A = pi +p 2 z 2 + 

P7qb 



(A14) 



A 
B 

C 

q 
t 



pi + Pit 

P6 + PlH 

Pn 



Pit 



l 



P3t 2 
P9t 



P5t, 



(l+P8tq) 2 
Viit 



Pwt, 



pnst, 



s(l+p 13 z 2 ) 3 
sin 2 (ip + P15), s = sin 2 (Lp+p 16 ), 
z = cos 2 (<p - ' 



Pi 



1 +P 3 z 



P5Z, 



B 



P6 



Piiqbt. 



6 1 



cos tp, 
In region II: 
A = pi 
B = P6 



P2Z 



Pi 



p 7 z 



1+P 3 z 

P8 
1 +P9Z 



■Pu)- 

■P5<P, 

■Pio<P, 



(A17) 
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Table 5. Coefficients pi for R p 



AC 



i 


I 


II 


III 


IV 




1 


0.049947 


-4.985248 


0.100241 


0.272905 




2 


-0.029006 


-0.025984 


0.005432 


0.058684 




3 


3872.363 


-0.007404 


-0.748377 


2.053694[- 


-3] 


4 


0.250385 


5.294455 


0.050631 


1.800867[- 


-7] 


5 


-0.245758 


-0.201654 


0.007900 


1.911708[- 


-8] 


6 


0.018241 


0.184431 


-0.032915 


0.052786 




7 


0.090256 


-0.139729 


-0.000768 


2.043824[- 


-5] 


8 


108.8302 


0.415562 


0.044312 


4.458912[- 


-8] 


9 


1.007326 


2.692073 


-0.697892 


1.101541[- 


-3] 


10 


0.061586 


-0.385832 


0.032534 


3.312811[- 


-4] 


11 


0.797695 


1.055347 


0.080109 


2.682799[- 


-4] 


12 


175.5965 


0.013667 


0.031994 






13 


9.306619 


-0.509106 


8.724039 






14 


-0.551550 


-0.267675 


2.982355 






15 


1.203014 


0.034585 


-0.062076 






16 


0.096598 










17 


-0.441039 











Table 6. Maximum values of v = u max in the fit ex- 
pressions and maximum relative error <5 max of the fits at 



Process 


I 


II 


III 




20 


23 


25 


<5max, % 


<10 


<15 


<20 


^AB : ^max 


15 


15 


25 


5max, % 


<20 


<20 


<50 


-^AC J ^max 


19 


21 


26 


<5max, % 


<15 


<26 


<20 


R p v 

AB ' max 


22 


13 


13 


5max, % 


<20 


<5 


<3 


"^AC ' ^ max 


23 


16 


13 


<5max, % 


<13 


<5 


<2 



C = pn +P12Z 2 
z = cos 2 <p. 
In region III: 
A = pi + p 2 (pz -+ 
B = p e + P7tpt + 

C = pn +Pl2<£ z 



P13 



1 +P14Z 



Pi 



1 +P3Z 

P&z 
1 + p 9 z 

Pl4 



Pl5f, 



P5<P, 
o rPlQif, 



t = sin 



1 +P13Z 
2irip \ 



0.321750554/ ' 



z = cos ip. 



Pwl+p 2 vl +P3vlv\ +P4V2, 



In region IV 
.4 
B 

C 



1 



P6 V 2 



-p 7 vl 



1 



p 9 v\+p w vl 



-rP5«f + 

■Puv%vf. 



(A18) 



(A19) 



(A20) 



Coefficients pi in (A17)-(A20) are given in Table 5. 
Calculation and fit errors do not exceed 5% for those val- 
ues for which i?^ c <; 10~ 6 . 

In Table 6 we give maximum values of v = v max of our 
fit expressions in regions I, II, and III, and maximum fit 
errors <5 max at v < w max in these regions. These maximum 
errors occur at v ~ w max , where the reduction factors are 
very small (and are thus unimportant for calculation of 
the total neutrino emissivity). At v > v max our fit ex- 
pressions are not reliable and we recommend to set the 
corresponding reduction factors equal to zero in computer 
codes. 



